Repository navigation
Transcriptome preset (-x tx), 2.3-2.8x faster than minimap2 for quantification - #27
Merged
Merged
Conversation
…hors dropped before chaining - -x tx: reads to a transcriptome for quantifiers (oarfish, NanoCount): minimap2 map-ont seeding (k15 w10, no occurrence range, -g5000) and up to 100 secondary alignments, with the dna preset's chaining and DP otherwise. - The ALT-contig count was recomputed for every read by scanning all targets (twice per read): 32 of 35 post-chain seconds on a 534k-transcript index. It is now counted once per index. - Anchors with no other anchor of the same reference/strand within the widest chaining distance are dropped before chaining (they can only form a one-anchor chain); the rescue's density test keeps counting all anchors. Output identical (DNA genome-wide sample 9bc7ff3e..., all golden hashes); DNA CPU -1.2%, transcriptome chaining -8%. - Region building copies chain anchors only for regions that survive the secondary filter. - Golden test tx_sam. HG002 200k reads vs minimap2 -ax map-ont -N 100 (GENCODE v49 transcripts): cDNA 185 -> 81 s, direct RNA 106 -> 37 s; same reads mapped; simulated reads at or above minimap2's accuracy; oarfish gene counts r = 0.9995 / 0.99997. Co-Authored-By: Claude Opus 5.5 <noreply@anthropic.com>
… transcriptome validation Co-Authored-By: Claude Opus 5.5 <noreply@anthropic.com>
Merged
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Add this suggestion to a batch that can be applied as a single commit.This suggestion is invalid because no changes were made to the code.Suggestions cannot be applied while the pull request is closed.Suggestions cannot be applied while viewing a subset of changes.Only one suggestion per line can be applied in a batch.Add this suggestion to a batch that can be applied as a single commit.Applying suggestions on deleted lines is not supported.You must change the existing code in this line in order to create a valid suggestion.Outdated suggestions cannot be applied.This suggestion has been applied or marked resolved.Suggestions cannot be applied from pending reviews.Suggestions cannot be applied on multi-line comments.Suggestions cannot be applied while the pull request is queued to merge.Suggestion cannot be applied right now. Please check back later.
Adds
-x txfor mapping direct RNA / cDNA to a transcriptome for quantifiers (oarfish, NanoCount): minimap2map-ontseeding (k15 w10), up to 100 secondaries, thednapreset's chaining/DP otherwise.HG002 200k reads vs
minimap2 -ax map-ont -N 100(GENCODE v49 transcripts): cDNA 185 → 81 s, direct RNA 106 → 37 s; identical number of unmapped reads; simulated reads land on their true transcript at least as often; oarfish gene counts r = 0.9995 / 0.99997. Report:docs/validation/2026-10-21-transcriptome/.General speedups (output identical: DNA genome-wide MD5
9bc7ff3e…, all golden hashes):Golden test
tx_sam; docs (README, usage, pipelines, benchmarks, changelog, validation index).🤖 Generated with Claude Code